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Abstract 

We update the study of the QCD corrections to direct J/ifi and T hadroproduction in association with a photon in the QCD-based approach 
of the Colour-Singlet (CS) Model. After comparison with the recent full next-to-leading-order (NLO) computation for this process, we 
provide an independent confirmation to the inclusive case that NLO QCD corrections to quarkonium-production processes whose LO 
exhibits a non-leading Pr behaviour can be reliably computed at mid and large by considering only the real emission contributions 
accompanied with a kinematical cut. In turn, we evaluate the leading part of the a'^a contributions, namely those coming from (J /if/, T) + y 
associated with two light partons. We find that they are dominant at mid and large Pr - This confirms our expectations from the leading Pj 
scaling of the new topologies appearing at NNLO. We obtain that the yield from the CS becomes one order of magnitude larger than the 
upper value of the potential Colour-Octet yield. The polarisation of the ^5 1 quarkonia produced in association with a photon is confirmed 
to be longitudinal at mid and large Pj . 
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1. Introduction 

For a long time, the many difficulties to correctly predict 
quarkonium-production rates at hadron colliders have been 
attributed to non-perturbative effects associated with chan- 
nels in which the heavy quark and antiquark are produced in 
a colour-octet state [1, 2, 3]. On the basis of leading-order 
(LO) (c^) calculations, it had been assumed that colour- 
singlet production channels give a small contribution at mid 
and large Pj-- The confusion most probably came from 
the fact that quantum-number conservation (J, P, C, and 
colour) prevents leading Pt scaling at LO and, in glue-glue 
production, at next- to-leading order (NLO) (a^)^. Thus, 
in contrast with the situation for many other observables, 
there is still a possibility of unexpectedly large (colour- 
singlet) contributions at large Pj at next-to-next-to-leading 
order (NNLO) (ff^ ). In addition, two unexpected features of 
fragmentation approximations have recently been revealed: 
heavy-quark fragmentation only dominates over the other 
topologies at very high Pt and gluon fragmentation may not 
dominate over double f-channel gluon exchanges. As a con- 
sequence, it seems that, when leading kinematic contribu- 
tions are correctly accounted for, colour-singlet production 
channels will play a more important role -if not the most- 
than in past analyses in the explanation of quarkonium- 
production observables and this suggests in turn that colour- 
octet fragmentation contributions may be less important 
than had previously been thought. 
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This became clear thanks to the several recent computa- 
tions of QCD corrections to quarkonium hadroproduction 
processes. The NLO corrections to the inclusive yield of 
J/t// and T were computed [4, 5] in the QCD-based ap- 
proach of the Colour-Singlet (CS) Model ^ [6]. Its polar- 
isation was in turn computed [8] at NLO. These computa- 
tions were recently complemented [9] by the addition of the 
real NNLO corrections -thereafter referred to as NNLO*-. 
It was then shown that there may be no need to incorpo- 
rate Colour-Octet (CO) transition (higher-v corrections of 
NRQCD) to describe the hadroproduction of T at the Teva- 
tron [10, 11, 12, 13]. In the case of the J/i// and i//' [14, 15], 
the CS contributions are significantly enhanced and brought 
very near the experimental data of CDF although the large 
Pt direct yield seems not to be fully accounted for in the 
case of the i^' for instance. As regards the CO channels 
for J/^ production, their NLO QCD corrections were re- 
cently computed in [16]. It was seen that they minimally 
affect both the Pt dependence and the normalisation of the 
partonic matrix elements at mid- and large-P^, thus also 
the value of the ^Sf^ CO Long-Distance Matrix Elements 
(LDME) fit to the data. Similarly, the polarisation pre- 
diction are not modified and remains in disagreement with 
the latest CDF measurements [17]. For recent reviews, the 
reader is guided to [1, 2, 3] along with some perspectives 
for the LHC [18] and a recent discussion on those afore- 
mentioned QCD correction computations [15]. 

Beside the studies of inclusive production, efforts are 
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Figure 1: Representative diagrams contributing to the hadroproduction of a J/ifr in association with a photon at orders a^a 
(a), ala (b,c), aia (d,e,f). See discussions in the text. 



being made to obtain improved theoretical predictions for 
complementary observables to the inclusive yield, such as 
hadroproduction of J/tfr and T (thereafter commonly named 
Q) in association with a photon [19, 20, 21, 22, 23, 24, 25, 
26]. Recently, the NLO corrections to the hadroproduction 
of (3 + 7 in the CS channel were computed by Li and Wang 
in [27]. As in the inclusive case, NLO corrections are sig- 
nificant in the large Pt region since new topologies appear 
with slower Pj falloff (see Fig. 1 (c)), in comparison to LO 
topologies Fig. 1 (a) and other NLO topologies as the loop 
corrections Fig. 1 (b). 

In [9], it was shown that some contributions to the 
inclusive yield coming from three-jet configurations, i.e. 
Q + jjj, such as those arising from gluon-fragmentation 
and "high-energy enhanced" (or double f-channel gluon ex- 
change) channels are important and clearly dominate over 
the other contributions at mid and large Pj. In the case 
which interests us here, (3 + y, those two channels appear 
at order a'^a - see for instance Fig. 1 (d) and (e,f)- and are 
expected to dominate over the Q + y yield at mid and large 
Pt. 

In this work, we therefore apply the same procedure as 
in [9] to evaluate their contributions. First, we will check 
that the NLO computation of Li and Wang [27] can indeed 
be reproduced by an evaluation of the real a^a contribu- 
tions (Q + J + one light partons) complemented by a cut- 
off on low invariant masses for any pairs of external light 
partons in the process. The procedure is explained in the 
next section. We will also check that the sensitivity of our 
computation on this cut-off dies away when Pj grows and 
that the theoretical uncertainty attached to its choice is typ- 
ically smaller than the ones attached to the choices of the 
renormalisation scale yu^ for as, the factorisation scale jj.f 
for the (collinear) parton distribution functions (PDF) and 
the heavy -quark mass nig. 

Having performed those checks, we will apply the same 
procedure for the evaluation of the contribution for Q + y + 
two light partons -namely the real NNLO contributions to 
Q + y production- arguing that they provide with a first reli- 
able estimate of the complete NNLO contributions to <3 + y 
production at large enough Pj- As regards the polarisation, 
the quarkonia directly produced in association with a pho- 
ton via those channels are mainly longitudinally polarised, 
as in the inclusive case. 



2. Cross section at NLO 

2.1. Inclusive case 

As we argued in [9], the NLO contributions to the in- 
clusive yield can be approximated at large enough Pj in 
a relatively simple and reliable manner by computing the 
ff^ contributions consisting in the production of a (3 with 2 
light partons (denoted j thereafter) on which we apply a cut- 
off on low invariant masses for any light parton pairs in the 
process. Computations of such cross sections can be done 
reliably using the automated generator of matrix elements 
MadOnia [28]. 

The underlying idea supporting this was twofold: 

• First, at large enough Pj, topologies with the leading 
Pj behaviour will dominate and those are wholly in- 
cluded in this subset of contributions (the produc- 
tion of a (3 with 2 light partons); 

• Second, this subset accounts for a physical process at 
Born level. Its contribution is therefore finite except 
for soft and collinear divergences. The purpose of 
the cut-off is to avoid such divergences by imposing 
a lower bound on the invariant-mass of any light par- 
ton pairs (sij). For the new channels (with a leading Pj 
scaling) opening up at , the dependence on this cut 
gets smaller for large Pj since no collinear or soft di- 
vergences can appear there. For other channels, whose 
Born contribution is at , the cut would produce log- 
arithms of Sij/s'^'". Those can be large. Nevertheless, 
they can be factorised over their corresponding Bom 
contribution, which scales at most as P^^. The sensi- 
tivity on s'™" is thus expected to vanish at large Pj. 

This insensitivity to the cut and the good agreement be- 
tween the NLO* (<3 + jj with a sij cut) and the full NLO 
result is recalled in Fig. 2 for the case of the inclusive J /if/ 
and Y(15) production. The gray band illustrates the sen- 
sitivity to the invariant-mass cut i™" between any pairs of 
light partons when it is varied from to 4m^ and 0.5m^ to 
2ml. ^'^ both NLO and NLO* computations, the value of 
all parameters were set to the same values. For the J/t//, 
we have m, = 1.5 GeV, \R(0)\^ = 0.810 GeV^ yu/ = 

ju, = = ^4m2 + Pi and BiiJ/if, fi+fx-) = 0.0588 
and, for the T(15), nib = 4.75 GeV, \R(0)\^ = 6.48 GeV^ 
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Figure 2: (a) full computation at NLO for J/ifr + X (dashed line) [4] vs. NLO* {J/i// + 2 light partons with a cut on Sij) (gray 
band) at = L96 TeV; (b) full computation at NLO for T(15) + X (dashed hne) [4] vs. NLO* (Y(15) + 2 light partons 
with a cut on s,y) (gray band) at = 1.96 TeV. See text for details, Br stands for the respective branching into dileptons. 
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Figure 3: (a) full computation at NLO for J/ip + y + X (dashed line) [27] vs. J/iJ/ + X + 1 light parton with a cut on i,y (gray 
band) at yfs ^ 14 TeV; (b) full computation at NLO for T(15 ) + y + X (dashed line) [27] vs. T(15 ) + y + 1 Ught parton with 
a cut on Sij (gray band) at = 14 TeV. The absolute value of the rapidity of both the Q and the y is limited to 3. 



= = yuo = y4OT^ + Pi and Br(T ^ yu+yu") = 0.0218. 
The parton distribution set used was CTEQ6JVI [29]. The 
yield becomes insensitive to the value of i™" as Pt in- 
creases, and it reproduces very accurately the differential 
cross section at NLO accuracy, both for the J/i// and T(15) 
case. 

2.2. Production in association with a photon 

In the more exclusive case Q + y, similar topologies are 
present with the same Pj scaling and we also expect to re- 
produce accurately the yield at NLO accuracy (a^ a) com- 
puted in [27] by computing the yield from the production 
of Q + y with one light parton with the invariant-mass cut 
i™" between any pairs of light partons, also referred to as 
NLO*. 



This is indeed the case. For instance, the differential cross 
section for Q-Fy at NLO accuracy from Li and Wang [27] is 
displayed in Fig. 3 and is very well reproduced by the NLO* 
computed for different values of 5™". When Pj grows, the 
latter becomes completely insensitive to the value chosen 
for i™". The same parameter values were used for Fig. 3 as 
for Fig. 2 and a was set to 1/137. 

This result is a clear and completely independent confir- 
mation of the validity of the reasoning initially given in [9] 
that NLO QCD corrections to quarkonium-production pro- 
cesses whose LO shows a non-leading Pj behaviour can be 
reliably computed at mid and large Pj by considering only 
the real emission contributions accompanied with a kine- 
matical cut. In turn, this reinforces our confidence that the 
impact of NNLO contributions can evaluated likewise by 



computing the NNLO* contributions, as done in the fol- 
lowing section. 

This also gives us confidence that much better Monte 
Carlo simulations of inclusive production at mid and large 
Pt could be achieved using NLO* and NNLO* partonic 
matrix elements, which can be interfaced [30] with event- 
generators such as PYTHIA [31]. In any case, they would 
give results much more reliable than simulations based on 
matrix elements for CS channels at LO only. 



3. Cross section and Polarisation at NNLO* 



3.1. Cross section 

Among the contributions appearing at a'^a, we find the 
topologies of Fig. 1 (d) (gluon fragmentation) and Fig. 1 
(e,f) ("high-energy enhanced" or double f-channel gluon ex- 
change), those exhibit new kinematical enhancements ap- 
pearing in higher-order QCD corrections. In other words, 
those provide us with new mechanisms to produce a high- 
er <3 with a -y with a lower kinematic suppression, still via 
CS transitions. They are therefore expected to dominate the 
differential cross section at NNLO accuracy in the region of 
large transverse momentum. 

Those are also entirely contained in the contributions to 
pp Q + y + jj, namely the real a'^a corrections, and we 
can follow the procedure validated in the previous section, 
by "simply" adding one light partons in the final state. 

The computation of pp Q+y+ jj at tree level is in prin- 
ciple systematic, but technically quite challenging: a dozen 
parton-level subprocesses contribute, most involving a few 
hundred Feynman diagrams. As done in [9], we follow the 
approach described in Ref. [28], which allows the automatic 
generation of both the subprocesses and the corresponding 
scattering amplitudes. 

The differential cross-sections for J/if) + y and T{lS) + y 
are shown in Fig. 4. The gray band (referred to as NLO*) 
corresponds to the sum of the LO and the real a contri- 
butions. The red (or dark) band (referred to as NNLO*) 
corresponds to the sum of the LO, the real a^^a and the 
real a'^a contributions. The a* a contributions in both 
case dominate over the yield at large Pj. The uncer- 
tainty bands are obtained from the combined variations 
O.SyUo < firj < 2po with for the J/tfr, nic - L5 + 0.1 
GeV and 1 < s'J'pil.S GtVf < 2 and, for the T(15), 
nib = 4.75 + 0.25 GeV and 0.5 < i™7(4.5 GeV)^ < 2. 

At the leading order in the heavy-quark velocity (v), the 
results for the radially excited states i/'(25'), T(25') and 
T(35) are readily obtained by changing |/?(a(0)p and the 
branching ratio into dileptons. 



3.2. Polarisation 

As regards the polarisation parameter A, it is computed 
by analysing the angular distribution {0) between the {'^ di- 
rection in the quarkonium rest frame and the quarkonium 
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Figure 4: Results for (up) Jjif/ + y and (down) T(15') + y 
from full NLO, NLO* and NNLO* contributions at = 
14 TeV. The theoretical-error bands for NLO* and NNLO* 
come from combining the uncertainties resulting from the 
choice of ^f, jj.^, and s™". The absolute value of the 
rapidity of both the Q and the y is limited to 3, > 3 GeV. 
A photon isolation cut (AR > 0. 1) was applied on the NLO* 
and NNLO* yields (see discussion in section 5.) 



direction in the laboratory frame. The normalised angular 
distribution /(cos 0) then reads 



/(cos 0) 



-(1 +^ cos^ 0), 



(1) 



2(A + 3) 

from which we can extract A bin by bin in Pj. 

Our results for the J/t// are shown in Fig. 5 along with the 
curves for the NLO*. Our predictions for the polarisation 
parameter A for the NLO* are in qualitative agreement with 
those of [27], i.e. the J/tfr's produced in association with 
a photon are dominantly longitudinal. We did not go fur- 
ther in the comparison since, for both numerical computa- 
tions, 5% precision in A could only be reached at a high cost 
of computing time and since the NNLO* is anyhow larger. 
As regards the latter, it confirms the trends of the NLO* 
(and thus NLO) results. This is not a surprise knowing the 
NNLO* results for the inclusive yield, differing essentially"* 



at least as far as the dominant channels are concerned. 
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Figure 5: Polarisation of the J lip produced in associa- 
tion with a photon at a/^ = 14 TeV up to the order of^or 
(NNLO*). Most of the uncertainties on A from the choice 
of and cancel. The uncertainty band of the NLO* and 
NNLO* result comes from the variation of the cutoff i™". 

in the replacement of the photon by a gluon. This replace- 
ment is indeed not expected to change results concerning 
the polarisations of the particles produced. 

4. Discussion of the results at NNLO* 

First, let us stress that although the uncertainty associated 
with the choice of the cut .s™" is somewhat larger than at 
NLO*, it is nevertheless smaller than the one attached to the 
mass, the renormalisation scale and the factorisation scale ^. 
The latter dependence is expected: on the one hand, we miss 
the virtual part at low Pj where it is sizable and where it is 
expected to reduce the renormalisation-scale dependence; 
on the other hand, the contributions dominating at large Pj 
are directly sensitive to the fourth power of as- Yet, the 
dependence is smaller than in the inclusive case [9] where 
five powers of as are involved. 

Second, we find that the subprocess gg — > Q+y+gg dom- 
inates, providing with more than two thirds of the whole 
yield in the case. In addition, we have checked that 
this fraction is slightly increasing with Pj and only weakly 
dependent on the value of the invariant mass cut-off of 
light partons i™", removing the coUinear and infrared diver- 
gences of gg ^ Q + y + gg. Another important observation 
is that the size of the yield from gg ^ Q + y + gg does not 
vary much when i™" is changed from 2.25 to 9 GeV^ in the 
y/i/f case for instance. This indicates that those divergences 
are not -after being cut- artificially responsible for a large 
part of the NNLO* yield. 

Third, it is likely that the largest part of this contribution 
is not from gluon fragmentation topologies, but rather from 



'Note that the uncertainty associated the choice of has a negligible 
impact on the final results compared to the other theoretical uncertainties. 



double f-channel gluon exchange ones*", keeping in mind 
that such a decomposition in terms of the corresponding 
Feynman graphs is not gauge invariant. There are a cou- 
ple of indications supporting this: 

• Such fragmentation contributions would be expected 
to provide transverse quarkonia if there were neither 
corrections due to the off-shellness of the fragment- 
ing gluons nor possible spin-flip contributions when 
radiated gluon energies (in Q rest frame) are not small 
compared to the heavy-quark mass. While a depolari- 
sation (in the helicity frame) is possible, it would be 
still far from the longitudinal polarisation computed 
here for Q produced via gg ^ Q + y + gg. 

• Processes such as qq' Q + y + qq' , which proceed 
uniquely via double f-channel gluon exchange, have 
the same Pt dependence as the process gg ^ Q + 
y + gg and the difference in normalisation is naturally 
accounted for by colour factors and the smaller value 
of the quark PDF compared to the gluon one at low 
X. Another similarity with gg — > Q + y + gg is that 
the polarisation of the yield from qq' Q + y + qq' 
is strongly longitudinal, for Pt larger than 5 GeV, as 
observed for the full NNLO* yield dominated by gg — > 
Q + 7 + 8g- 

The previous discussion can in fact be extended to the 
inclusive case studied in [9]. Indeed, the results then ob- 
tained at NNLO* appeared to be significantly higher than 
the ones using the fragmentation approximation. Further, 
the polarisation of the yield from gg ^ Q + ggg processes 
happened to become more and more longitudinal at large 
Pt in contradiction with the expected polarisation from a 
fragmentation channel. This lets us think that, for both pro- 
cesses Q + y and Q + X at NNLO*, the yield may mainly 
come from double f-gluon channel exchanges appearing for 
the first time at this order. 

A careful kinematical analysis of the yield from gg — > 
Q + y + gg would certainly be helpful. Yet, it would be 
highly computer-time demanding, especially to obtain a dis- 
tribution^ of the relative momentum between the J/tfr and 
its closest gluon precise enough to unequivocally attribute 
the most part of the yield to the double f-channel gluon ex- 
change channels rather than to the fragmentation ones. In- 
deed, it concerns the most complicated process with a cou- 
ple of hundreds of diagrams. This work is left for a fu- 
ture analysis and would certainly be expediently done along 
with the one for the inclusive case. 



In the inclusive case, those have been previously discussed [32] in the 
k, factorisation formalism. See [33] for a recent application to T hadropro- 
duction. 

^The evolution of the distribution for different Pt would be even better. 



5. Phenomenology 

5.1. Photon detectability 

In order to detect the photon, we evidently have to im- 
pose that it possesses a finite transverse momentum for it 
not to go in the beam pipe and that it is isolated to avoid 
misidentifications with Bremsstrahlung radiations. Yet, iso- 
lation criteria highly depend on the detector potentialities. 
The determination of optimum values for and for an 

isolation criterion for the y improving the signal over back- 
ground ratio is beyond the scope of this theoretical analy- 
sis. We have therefore decided to apply minimal values for 
those constraints, which would be further increased after a 
full detector simulation. Along the same lines, it is worth 
noting here that isolation cut for CO mediated signal (a pri- 
ori suppressed, as discussed below) can only be imposed in 
such simulation. Indeed, such a cut requires a simulation of 
the hadronic activity following the discolouration of the CO 
pairs as done e.g. in [34]. 

At LO, the Pj requirement is trivially satisfied for a Q 
with a finite P^ since is balancing P^. As discussed by 
Li and Wang [27], this is not automatically the case at NLO 
and a minimum cut has to be applied and it aff'ects the 
yield vs P^ up to roughly three times the value of this cut. 
Typically the cut P^„„„ = 3 GeV applied here affects the 
yield up to = 7 - 10 GeV. As regards the isolation cut, 
we excluded any events with partons (q,g) within a cone 
A/? - ^[^rf^¥~K(j? - 0.1 from the photon [35]^. As noted 
above, those are minimal cuts. 

In our case an analysis of different P^-cut values is less 
relevant than in [27] for two reasons: first, our computa- 
tion cannot be reliably extended to low Pj and, second, the 
main contributions -the double f-channel gluon exchange 
and gluon-fragmentation channels- create a y with a sim- 
ilar momentum as that of the Q (see the graphs of Fig. 1 
(d) & (e)). Overall, such a kinematical cut insuring its de- 
tectabiUty does not affect the NNLO* yield for Pj larger 
than 10 GeV. 

Beside the problem of the photon detectability, we have 
to make sure that other processes will not contribute to 
the yield of Q + y. In the Jjili case, we expect the non- 
prompt background to be properly subtracted (by vertex- 
displacement method for instance). However, in this case, 
the ratio "non-prompt over prompt" at very large Pj is ex- 
pected to be lower than in the inclusive case since the pho- 
ton emission is required for both processes, whereas in the 
inclusive case the B feed-down can proceed at a lower order 
in as ■ Similarly, the xq feed-down will not be significant 
(except in the region where the invariant mass of the pair 
<2 + y is close to the mass of the xq) since suppressed by 



On the way, note for completeness that this cut avoids the QED singu- 
larities that may appear in the processes qg — > qgyQ when the y is emit- 
ted by the external quark. Those are anyhow sub-dominant topologies of 
suppressed quark-gluon initiated contributions and could have safely been 
neglected. 



and by the branching while having the same suppres- 
sion than for the Jjij/ for similar topologies. As regards the 
feed-down of the radially excited states ^S\, they are read- 
ily accounted for by constant multiplicative factors to the 
cross section of the state fed in : ~ 1 .4 for the i^(25 ) into 
J I if/, -1.1 for both the T(25) into T(15) and T(35) into 
T(25) ^ 

5.2. Colour-octet yield 

Now, let us discuss the possible colour-octet contribu- 
tions, ignoring first the modifications, induced by the QCD 
corrections to CO and CS contributions, of the CO LDME 
extracted from the Tevatron data. As discussed in [27], a 
quick comparison with the results obtained in [26] shows 
that the NLO CS and the LO CO yields are of the same 
order at the LHC, with roughly the same Pj dependence. 

Yet, contrary to what is claimed in [27], higher-QCD cor- 
rections are expected to impact more on the CO partonic 
matrix elements mQ + y than in the inclusive case [16], for 
which NLO corrections seem unimportant at large Pj- In- 
deed, the difference between the two cases lies in the fact 
that in Q + y Pj!^ topologies, such as fragmentation ones, 
are not opened at the order aa^ as studied in [26] (only 
one initiated by quark, and thus sub-dominant at low x, 
qq ^ y + {QQ)^^\ with a quark in the f-channel shows a 
Pjf* scaling). 

In the light of the NNLO* results, we know that when 
the real a^a contributions are taken into account, the CS 
contributions will be increased by about one order of mag- 
nitude compared to the CS at NLO. To be sure that CS con- 
tributions are then dominant over CO ones at higher orders, 
we have to check that fragmentation contributions from CO 
opening at aa^^ are small. 

This is indeed the case since they are mediated by ei- 
ther or 3p{f^ C = +1 CO states and their coiTe- 
sponding LDMEs are know to be severely constrained when 
higher-QCD corrections are taken into account for the anal- 
ysis of the inclusive data at low- and mid-Pj- at the Teva- 
tron [16, 36], at various fixed-target experiments [37] and 
more recently at the B factories [38, 39]. Finally, a fur- 
ther suppression is expected from the increase of the CS 
yield from NLO and NNLO* contributions [9]. Using an 
upper conservative value of the LDME combination M^** = 
<0|6)^/'^[i5{,^^]|0>H-4<0|O-'/^[3p{fl]|0) = IxlO-^GeV^ (k = 
3.5), we analysed the two extreme cases «0|(9^/^[3p{fl]|0) = 
(case I) and {0\O^''''[^S^q^]\0) = (case II)) to evaluate 
those CO aa^ contributions via pp ^ j + (cc)** *-"^' + y 
using the same parameter values as the central NLO curve 
and i™" - ni^c- Pt - 10 GeV, the case I gives a con- 
tribution about 5 times smaller than the central aa^^ (NLO) 
CS values while the case II is further below'" by another 



'The other feed-down factors can be neglected in view of the low accu- 
racy at which we know |Rq(0)P and the respective branchings from which 
they are derived. 

'"This indicates on the way that this process is sensitive to a different 
linear combination of C + 1 CO LDME than the inclusive production. 



factor of 5. At - 40 GeV, case I is less than a factor of 2 
larger than the central NLO values, thus still clearly below 
the NNLO* yield, while case II is still a bit suppressed and 
comparable to the central NLO values. A dedicated analysis 
including a global fit of the C = +1 CO LDMEs and taking 
into account DGLAP radiation of the fragmenting gluon is 
beyond the scope of this work but would confirm that aa^ 
CO corrections to Q + y matter only at large Pj and are 
most likely not large enough to challenge the dominance of 
CS contributions to this process. 

This dominance is expected since the photon has to be 
emitted by quarks. If the hard scattering part is gg gg, 
the photon can only be emitted by the heavy quarks with 
topologies similar to CS channels. At a^a, the gluon-fusion 
CO channels do not scale like Pj"^. One has to go to of^of 
to have the first gluon fragmentation channels initiated by 
gluon fusion. In this case, as discussed above, the gluon 
fragments into a photon and a CO C = + 1 (whose LDMEs 
are severely constrained) ' ' . Last but not least, contrary to 
the double f-channel gluon exchanges dominating the CS 
yield, the latter fragmentation CO topologies will system- 
atically produce nearly collinear y and hadrons from the 
discolouration of the CO state. Those events are therefore 
likely to be rejected by experimental photon-isolation cuts. 

Finally, while contributions from i-channel cut to J/tfr 
production could appear in the low Pj- region as in the inclu- 
sive case [40] (the final state gluon of the gg — > J/t//g being 
simply replaced by a photon), in the large Pj region, the 
colour-transfer-enhancement mechanism discussed in [41] 
is not expected to matter for the present process. In the in- 
clusive case, gg ^ J/iff + cc may account for a significant 
part of the yield at very large Pj. This indicates that topolo- 
gies for which colour transfer could occur are not much 
suppressed and those could impact on the inclusive yield. 
In the present case, we have checked that the contribution 
to J/i// + y from the partonic process gg — > Jjif/ + cc + y 
is sub-dominant and found that it is one order of magni- 
tude smaller than the process from gg — > J/iff + y + gg 
with the same Pt dependence between 20 and 50 GeV. It 
is therefore quite unlikely that colour transfers, acting on 
topologies with at least 3 heavy quarks, be visible inJ/i/z+y 
yields. 

6. Conclusion 

In conclusion, we have computed the real next-to-next- 
to-leading order QCD contributions to the hadroproduction 
of a y/i/' -H y and T + y via colour singlet transitions along 
the same lines as [9] for the inclusive case and argued that 



' ' If one required that exactly balance in order to select the LO 
CSM, aiming at the extraction of the gluon distributions for instance, one 
would have to refine the analysis to take into account NNLO* corrections 
in the CS channels and the CO yield [25] which would then be dominated 
by gg^ g* ^ (22)"*' + 7 via 'S™ or V™ if the C = +1 CO channels 
are not too much suppressed. 



it provided a first reliable evaluation of the corresponding 
yield at NNLO accuracy. 

Indeed, prior to that, we have shown that the full NLO 
evaluation of Li and Wang [27] is very accurately repro- 
duced for Pt > 5 GeV by the sole evaluation of the real 
emission contributions up to order a^^a, namely the NLO*. 
In the inclusive case [9], a similar observation was made and 
motivated the study of the NNLO* contributions to evaluate 
the yield at NNLO accuracy. We have thus reach a com- 
pletely independent, but similar, conclusion for the produc- 
tion of a quarkonium in association with a photon. 

A priori, integrating the amplitudes for such real emis- 
sion contributions leads to divergences in some phase-space 
regions. Those can be avoided by imposing a minimal in- 
variant mass between the external light partons, while in a 
complete calculation those divergences are canceled by the 
virtual corrections. Yet, at NLO, the latter scale as Pj^ and, 
at NNLO, as Pj^, hence are suppressed compared to respec- 
tively the real NLO part scaling as Pj^ and the real NNLO 
part scaling as Pj'^. This explains why we can neglect the 
virtual corrections at mid and large Pt and use a cut-off 
on which the results become insensitive when Pt grows, 
and also why the sole NLO* reproduces the full NLO accu- 
rately. 

We have also shown that the diff'erential cross section in 
Pt at NNLO* is one order of magnitude larger than at NLO 
for large Pt -as expected from their Pt scaling-. We have 
also identified the process responsible for the most part of 
it, i.e. gg ^ Q + y + gg. Although, the distinction between 
the gluon-fragmentation and the double f-channel gluon ex- 
change graphs cannot be carried out in gauge invariant way, 
we provided some hints that the second type of topology is 
dominating, similar to the inclusive case [9]. One of this 
hint is the polarisation of the quarkonia produced in associ- 
ation with a photon. 

Indeed, we have computed the polarisation parameter A 
of the NNLO* yield which is negative, indicating a longitu- 
dinally polarised yield. This confirms the trend observed at 
NLO and hints at the dominance of double f-channel gluon 
exchange contributions. 

When the NNLO* contributions are incorporated in the 
CS yield, it becomes one order of magnitude larger than the 
potential CO yield, which would mainly produce transver- 
sally polarised Q. The measurements of the cross section 
for the production of Q + y would directly measure the size 
of the CS without being sensitive to the non-perturbative 
CO parameters. This is a complementary case to the study 
of J/iJ/ + cc and T + bb as discussed in [5, 42, 15]. 

Similarly to the analysis of the inclusive case [9], this 
analysis cannot be extended to too low Pt, where the ap- 
proximations on which it is based no longer hold. One way 
to improve the predictions could be achieved by merging 
the matrix elements with parton showers using one of the 
approaches available in the literature [43], or by performing 
an analytic resummation [44] . 

Finally, the results presented here strongly support the 
procedure used for the very first evaluation of the inclu- 



sive yield at NNLO accuracy [9] and which showed an 
agreement with the Tevatron measurements. This in turn 
confirms that much better Monte Carlo simulations than 
the ones based on matrix elements for CS channels at LO 
only are now possible at mid and large Pj. This could 
be achieved using NLO* and NNLO* partonic matrix ele- 
ments generated by MadOnia [28, 30], for the process stud- 
ied here pp — > Q + j + X, but also for the inclusive mea- 
surements to be performed at the LHC. 
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